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ABSTRACT 

We present Wisconsin Ha Mapper observations of ionized gas in the Smith Cloud, a high velocity 
cloud which Lockman et al. have recently suggested is interacting with the Galactic disk. Our Ha 
map shows the brightest Ha emission, 0.43 ± 0.04 R, coincident with the brightest H I, while slightly 
fainter Ha emission (0.25 ± 0.02 R) is observed in a region with H I intensities < 0.1 times as bright 
as the brightest H I. We derive an ionized mass of > 3 x 10 6 M Q , comparable to the H I mass, with 
the H + mass spread over a considerably larger area than the H I. An estimated Galactic extinction 
correction could adjust these values upwards by 40%. Ha and [S II] line widths towards the region of 
brightest emission constrain the electron temperature of the gas to be between 8000 K and 23000 K. 
A detection of [N II] A6583 in the same direction with a line ratio [N II]/Ha = 0.32 ± 0.05 constrains 
the metallicity of the cloud: for typical photoionization temperatures of 8000 — 12000 K, the nitrogen 
abundance is 0.15 — 0.44 times solar. These results lend further support to the claim that the Smith 
Cloud is new material accreting onto the Galaxy. 

Subject headings: Galaxy: evolution — Galaxy: halo — ISM: abundances — ISM: clouds — ISM: 
individual (Smith Cloud) 



1. INTRODUCTION 

The infall of gas is a significant driver of the evolution 
of galaxies. Without infall, there is insufficient gas in the 
interstellar medium of the Galaxy to sustain the observed 
star formation rate for more than ~ 1 Gyr. Also, the 
narrow distribution of mctallicities of long-lived stars in 
the solar neighborhood implies that the metallicity of 
the ISM is roug hly constant over tim e (the so-called "G- 
dwarf problem" ; |van den Bergh 1962| . In the Milky Way, 
high velocity clouds (HVCs) provide direc t observational 
evidence for infalling low metallicit y gas (Wakker & van 
Woerden|1997| IWakker et al.|2008|» 
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i'he Smith Cloild" is an excellent example of the active 
infall of material. Also called the Galactic center posi- 
tive (GCP) complex, the Smith Cloud is an HVC with a 
radial velocity near +100 km s -1 with respect to the lo- 
cal standard of rest ( |Smith|1963 ). Lockman et al. (2008 



(Wakker et al. 2008). Therefore, the position and trajec- 
tory of the cloud are both well known: the head is moving 
towards the plane at 73 ± 26 km s" 1 , and the portion of 
the cloud mapped by L08 has an H I mass of > 1O 6 M0. 

In this paper, we present Wisconsin Ha Mapper 
(WHAM) observations of ionized gas in the cloud, includ- 
ing the first spectroscopic Ha map and pointed, velocity- 
resolved, high signal-to-noise-ratio spectra of the Ha, 
[N II] A6583, [S II] A6716, and [O III] A5007 lines. 

2. OBSERVATIONS 

WHAM is a dual-etalon Fabry-Perot spectrometer cou- 
pled to a siderostat, designed to observe very faint opti- 
cal emission lines from diffuse gas. The instrument (de- 



gas. 

2003[ ) integrates all emission 
field of view, sacrificing spatial information 



scribed by |Haffner et al. 
within its 1 

to obtain spectra with a resolution of 12 km s _1 over a 
hereafter L08) presented a survey of H I 21 cm emis- 200 km s _1 -wide window with a sensitivity below 0.1 R in 



sion from the cloud using the Gr een B ank Telescope, 
and Bland-Hawthorn et al. (19981 and Putman et al. 
( 2003 hereafter P03) identified weak Ha and [JN II] emis- 
sion associated with the cloud. 1 Like many HVCs, the 
cloud has a cometary morphology with a bright tip and 
a diffuse, trailing tail. L08 showed that the cloud is in- 
teracting with the Galactic disk and used this informa- 
tion to obtain a kinematic distance of 11.1 — 13.7 kpc, 
consistent with a stellar absorption distance constraint 



a 30 s exposure. 2 We calibrated the geocentric velocity 
of the Ha emission with a fit of the geocoronal Ha line; 
this calibration is accurate to ~ 1km s . Other emis- 
sion lines are calibrated based on the instrument configu- 
ration and are accurate to within a few km s _1 (Madsen 
2004|. The observations presented here were obtained 
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'Throughout this paper, we use "H I" to refer to both neutral, 
atomic hydrogen gas and the 21 cm emission from the gas. We use, 
e. g., "H+" and "N+" to refer to other ions (ionized hydrogen and 
singly ionized nitrogen, respectively), with "Ho" referring to the 
Balmer-c? emission from recombining ionized hydrogen, H + . Other 
ionic species notations refer to the emission lines from the gas (e. 
g. "[N II]" refers to forbidden emission from N+ gas). The H I 
21 cm intensity is proportional to the H I column density, whereas 
the Ha intensity is proportional to the line integral of the square 
of the H + density. 



with WHAM while it was located on Kitt Peak in Ari- 
zona. 

2.1. Pointed observations 

In 2007 October, we obtained spectra towards the tip 
of the cloud using the following sequence. 120 s "ON" 
observations of the tip (I = 38.6°, b = —13.1°) were al- 
ternated with 120 s "OFF" observ ations (chosen based 



on the Bland-Hawthorn et al. 1998 H I map) for a total of 
360 s on source for Ha and 480 s on source each for [N II] , 
[S II], and [O III] and equal times off source. WHAM 
spectra have a substantial contribution from the bright 

2 1 Rayleigh = 10 6 /47r photons cm -2 sr _1 s . 
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TABLE 1 

Gaussian fit parameters of deep spectra 



50 100 

LSR velocity (km/s) 

Fig. 1. — Deep spectra of the the Smith Cloud tip (the blue X 
Fig. [2|. Data are shown as diamonds with error bars; the large 
uncertainties in the Ha spectrum near «lsr = —20 km s — 1 are due 
to the subtraction of the relatively bright geocoronal emission. A 
multiple-component Gaussian fit to the data is shown as a solid 
line with each component plotted as a dashed line; fit parameters 
are listed in Table m The average of L08 H I 21 cm spectra over 
the WHAM beam is shown in the bottom panel. 

geocoronal Ha line and numerous fai nter, mostly uniden- 



tified atmospheric emission lines. Hausen et al. (20021 
constructed an empirical template of the faint lines near 
Ha, which are stable relative to each other but vary in 
intensity, partly dependent upon zenith angle. Using this 
technique for other lines, we fit the OFF spectra of each 
line to a single atmospheric template, which we then sub- 
tracted from the ON spectra to obtain spectra consisting 
purely of astronomical emission. Summed, atmosphere- 
subtracted ON spectra for each line are shown in Fig- 
ure Q] 

We fit the atmosphere-subtracted spectra with three 
Gaussians convolved with the instrument profile: one 
for local gas near Wlsr = +10 km s _1 , one near i>lsr = 
+40 km s _1 , corresponding to the Sagittarius Arm of the 
Galaxy, and one for the Smith Cloud near +100 km s _1 . 
Solutions are listed in Table [l] with our adopted line 
widths and intensities listed in the bottom of the table. 
We first performed an unconstrained fit. Blending of the 
emission from the Smith Cloud and the relatively wide 
Sagittarius Arm component is the dominant source of 
uncertainty in constraining the multiple-Gaussian fits, 
although the width and intensity of the Smith Cloud 
component were well-defined in all three lines. The rel- 
atively narrow [S II] line provides the best constraint on 
the mean velocity of each component, so we fixed the 
Smith Cloud velocity in the Ha and [N II] spectra to the 



Mean (km s 1 ) 


FWHM (km s- 1 ) 


J(R) 




Ha, Smith mean fixed to [S II]. 


Xr = I- 2 - 


+11.2 ±0.5 
+43.2 + 4.5 
+96.2 a 


30.3 + 2.3 

58.3 + 7.5 

33.4 + 2.4 


1.06 ± 0.20 
1.21 ± 0.22 
n + n nzL 


[N II], Smith 


mean fixed to [S II]. 


Xl = 0-5. 


+13.3+ 1.0 
+40.2 + 5.7 
+96.2 a 


26.9 + 3.2 
48.2 ± 7.5 
24.4 + 3.7 


0.54 + 0.17 
0.71 + 0.18 
0.14 + 0.02 


[S II] 3-component free fit. Xr 


= 0.6. 


+11.3 + 0.8 
+43.5 + 4.0 
+96.2+ 1.1 


22.9+ 1.6 
39.1 ± 8.2 
19.7 + 3.3 


0.62 + 0.07 
0.35 + 0.07 
0.14 + 0.02 


[O III], Smith 


mean fixed to [S II] 


X 2 r = 0.3. 


+13.3 a 
+43.5 a 
+96.2 a 


25.0 a 
45.0 a 
25.0 a 


b 
b 

0.026 + 0.014 


Adopted values for the Smith Cloud tip 


Ha 
[N II] 
[S II] 

(o mi 


33.4 + 2.4 
24.4 + 3.7 
19.7 + 3.3 

25 a 


0.43 ± 0.04 c 
0.14±0.02 c 
0.14±0.02 c 
< 0.07 (3a) 



a Values with no uncertainty were fixed in the 
fit. b [O III] components at i>lsr < +50kms — 1 ac- 
count for excess residual emission in the ON— OFF 
spectrum and have no physical significance. Intensi- 
t ies a re not corrected for extinction, as discussed in 

§m 

+96.2 km s _1 of the [S II] line and performed a second 
fit. The uncertainties are statistical and do not account 
for systematic effects; for example, any unidentified com- 
ponents in the spectrum would cause the widths in our 
fit to be overestimated. 

For the [O III] spectrum, we have no suitable atmo- 
spheric template for the weak emission at the Smith 
Cloud velocities. Figure [T] shows an ON— OFF spectrum 
for this line; no atmospheric template is subtracted. We 
fit the spectrum by fixing the velocities of all three com- 
ponents to those observed for [S II] and the widths at 
values similar to the observed width of the [N II] com- 
ponents. This yields a best fit [O III] intensity in the 
Smith Cloud tip of 0.026 ± 0.014 R, with the reduced 
X 2 parameter Xr = 0.34. However, fitting a flat line to 
the spectrum (no [O III] emission) yields nearly as good 
a fit, with Xr = 0-40. This result is a 3tr upper limit 
of 0.07 R, fainter than a previous [O III] upper limit of 
0.12 R ( |Bland-Hawthorn et al.||1998[ ). 



2.2. Maps 

In 2007 May and June, we obtained a map of the Ha 
emission from the Smith Cloud r egion using WHAM with 
the "block" mapping technique (| Haffner et al.|2003| . Ob- 
servations of each sightlinc within a 7" x 7 block were 
obtained in succession with 60 s exposures for each point- 
ing. One-degree beams were observed on a grid with 
Al = 0.98°/ cos b and A6 = 0.85°, undersampling the 
image. We used an average spectrum of all pointings 
within the block to characterize the level of the atmo- 
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Fig. 2. — Map of the Smith Cloud in Ha in Galactic coordinates. Intensities are derived from the area of a Gaussian fitting the Ha 
spectrum centered at +96.2 km s _1 with a width of 30 km s . Greyscale contours of H I column density in the L08 data, integrated from 
+70 to +145 km s _1 , are shown at 2 (dark grey), 10 (light grey), and 20 (white) xlO 19 cm -2 . Yellow contours show intermediate velocity 
H I integrated from +50 to +60 km s — 1 at 4 X 10 19 cm -2 . A blue X sign indicates the postion of the on-source deep exposures, green X 
signs show the off-source positions, and green + signs show the positio ns of the P03 spectra . The bright Ha emission near (35°, —6°) is 
the edge of a low-extinction window into the inner Galaxy, discussed by Madscn & Reynolds (2005). The scale bar assumes D = 12.4 kpc. 
The green line is discussed in § [4] 

the pointed observation. 

In each fit, the means of the lower- velocity compo- 
nents were allowed to vary freely. Not all beams in the 
map contain emission at the velocity of the Smith Cloud. 
We tried fitting the full profile allowing the mean of the 
+96.2 km s _1 component to vary, but in sightlines with 
low or zero emission from the Smith Cloud, the lowest x 2 
fit occurred with that component fitting a residual from 
lower-velocity Galactic emission instead. To recover the 
most accurate intensity distribuion of the ionized com- 
ponent of the Smith Cloud, we fixed this component to a 
mean of +96.2 km s _1 to trace even the faintest Ha emis- 
sion associated with the cloud near this velocity. The 
narrower [S II] line is better suited to a full exploration 
of the velocity field in the Smith Cloud, but our current 
data set only has [S II] data along one sightlinc. 

2.3. Extinction 

We estimate the extinction towards the cloud using 
the mean H I column density in the direction of the 
WHAM observations integrated over —20 to +50 km s _1 , 
(iVjj j) = 9 x 10 20 cm' 2 ; based on its velocity, this emis- 
sion should be predominantly foreground. Assuming the 
mean ratio of color excess to H I column density found by 
Bohlin et al] ( |1978[ ), we estimate E{B -V) = 0.15 mag. 
If the dust fol lows t he extinction curves determined by 
( |1989[) with the standard diffuse ISM value 
myyEJB -V) = 3.1, this yields a total ex- 
tinction at Ha of A(Ha) = 0.38 mag. The scatter in the 




Fig. 3. — Smoothed map of the Smith Cloud in Ha, a portion 
of the area shown in Fig. [2] Contours identify the chosen H I 
boundaries of the clouds identified in Table [2] 

spheric template, which we then subtracted from each 
spectrum. We fit each atmosphere-subtracted spectrum 
with the sum of (1) a single Gaussian to account for the 
geocoronal Ha and (2) three Gaussians near LSR veloc- 
ities of +10, +40, and +96.2 km s _1 , as observed in the 
deeper, pointed observations. The widths of all three 
components were fixed to 30 km s _1 , typical of interstel- 
lar Ha emission and similar to the derived line widths in 



Cardelli et al 



foreground H I column density leads to a standard de- 



4 



viation of (JAVRa.) — 0.11 mag. This estimate is similar 
to that of P03 using COBE/DTRBF, dust emission data, 
resulting in an increase in the Ha intensity of about 40%. 
However, because the comparison of H I column densi- 
ties and extinction is uncertain and the line ratios of Ha, 
[N II], and [S II] are not affected by extinction due to the 
similar wavelengths of the lines, we do not attempt to ap- 
ply an extinction correction to our data. In the future, 
we will obtain H/3 spectra of the cloud to more accurate ly 
measure the Ha extinction (Madscn & Reynolds 2005). 



3. RESULTS 

In the Ha, [Nil], and [S II] pointed spectra in Fig- 
ure [T] the Smith Cloud is clearly detected near Ulsr = 
+100 km s _1 . Emission from the warm ionized medium 
(WIM) in the Saggitarius Arm is evident near vlsr — 
+40 km s _1 . In Figure [2] we show a map of the Ha 
emission from the Smith Cloud. For each spectrum, we 
derived intensities from the Gaussian fit described above. 
The region of brightest H I emission — the "tip" identified 
by L08 — is also brightest in Ha, with the pointed obser- 
vations indicating an intensity of 0.43 ± 0.04 R. How- 
ever, a separate component of the cloud near (l,b) = 
(38°, -19°) (labeled "clump A" in Fig. [3) that is present 
but faint in H I has an Ha intensity of 0.25 ± 0.02 R. 

Along a line from roughly (38°, -11°) to (43°, -15°), 
the +100 km s -1 H I and Ha emission both cut off 
against a narrow ridge; outside this ridge is bright 
intermediate- velocity H I emission, identified by yellow 
contours in Figure [2] L08 suggest that this is likely 
material that has been ram-pressure-stripped from the 
cloud. At the latitude of the Smith Cloud, the sight- 
lines pass through the Satittarius Arm at z ~ — 1 kpc 
and «lsr ~ +40 km s _1 . Due to the 1.0 — 1.8 kpc scale 
height of the warm ionized medium (e. g 



Haffner et al. 



1999| |Gaensler et al ||2008| |Savage fe Wakker||2009p and 
the < CT4kpc s cale height of the warm neutral medium 
(Ferriere 20011, any Ha emission associated with the 
intermediate- velocity H I emission cannot be separated 
from the warm ionized medium Ha emission from the 
foreground Sagittarius Arm at similar velocities. 

There is considerable pixel-to-pixel variation in the 
Ha emission on the 1° scales probed by WHAM. P03 
reported a 20% variation in Ha intensity between two 
5' beams spaced 0.5° apart, and the L08 21 cm data 
show structure on scales comparable to their 10' resolu- 
tion. Therefore, unresolved Ha emission structure may 
be present in the data presented here. 

3.1. Mass 



IWakker et al. 



Smith Cloud to 



( 2008 1 estimated the H + mass of the 
loeJUTI- 1.9) x 1O 6 M based on the Ha 
detections in two directions reported by P03. With the 
cloud mapped in Ha, we refine the mass estimate using 
one of two assumptions: (1) the H + is concentrated in 
an ionized skin, or (2) the H + is fully mixed with the H I 
gas. For each case, we derive the emission measure from 
the Ha intensity (EM = fn 2 e ds = 2.75T 4 °- 9 / Hq pccm" 6 , 
where T4 = T(10 4 K)~ 1 and In a is t he Ha intensity 
expressed in Rayleighs; |Reynoldsj|1991| ), assuming T 



10 K (see § 3.2 1. If a uniform Lyman continuum flux 



incident on the cloud maintains the ionization and the 
cloud is optically thick to Lyman continuum photons, 



the emission measure is constant across the cloud as seen 
from the source of the ionizing radiation. 

We es timate th e mas s treating several clumps sepa- 
rately (§ |3.1.1 and 3.1.2 1 and considering the entire cloud 
as a contiguous entity (§ 3.1.31. We first identified each 



clump, shown in Figure [3J based on localized Ha emis- 
sion and assigned WHAM beams to the clump by eye. In 
each case, there is a region of corresponding enhanced H I 
emission, so we then chose a minimum H I column den- 
sity, Min iVjj j, for each clump which defines the region 
of enhanced H I. These values are noted in Table [2) The 
"tip" — which contains the pointed observations discussed 
above — and "tail" are in regions of relatively bright H I 
emission, with Min iVjj j ~ 10 20 cm" 2 and are along the 
major axis of the cloud. "Clump A" and "clump B" are 
in regions with fainter H I, Min iVjj j = 2 x 10 19 cm~ 2 
and are each separated from the major axis of the cloud 
by several degrees. These four clumps were chosen sub- 
jectively based on their enhancement relative to the dif- 
fuse emission. Among these clumps, the tail is most sim- 
ilar in intensity (in both Ha and H I) to the nearby 
surrounding emission from the cloud. 

3.1.1. Ionized skin 

First, we assume that each clump consists of uniform- 
density neutral hydrogen with a fully ionized skin of the 
same temperature and pressure. In this case, the electron 
density in the skin n e — n a /2, where n is the density 
in the neutral clump. This condition is established when 
the ionized gas within the skin has had time to come 
to pressure equilibrium with the neutral clump. Alter- 
natively, if this gas has not had time to reach pressure 
equilibrium, the density in the skin could be as large as 
that in the neutral clump (n e = no). We derive physical 
parameters for an ionized skin model with these two lim- 
iting cases below. We further assume that the depth of 
the neutral cloud along the line of sight is comparable to 
the projected width, Ljj j. In the skin, the path length is 
derived from one of these density assumptions combined 
with our Ha observations: Lh+ — EMn" 2 . 

The mass of the ionized gas in the cloud within a solid 
angle f2 is lAm,nn e D 2 HL i i+ , where D is the distance 
to the cloud. The mass of a hydrogen atom is mn; the 
factor of 1.4 accounts for helium. The ionized mass in 
the cloud within one WHAM beam is thus 



Ah 



H + 



1.27 x 10 J 



D 



EM 



12.4 kpc J peem 



V cm 6 J 



For each clump, we derive the neutral gas density as 
(no) = (A+j i)/ijj j, where (ATjj t) is the mean H I col- 
umn density over the region where Ajj j > Min A+r-r. 
We then use n e and the observed EM with equation |l]) 
to calculate the H + mass in each beam. The H + mass m 
each clump is then the sum of the masses in the beams 
identified with the clump. The resulting total H + masses 
and path lengths are listed in Table [2] for both the equal- 
temperature-and-pressure and the equal-density cases. 
Using this ionized skin assumption, there is w 1 x 10 5 Af Q 
of ionized gas in the four clumps. Note that the tip and 
tail, which contain the brighest Ha and 21 cm emission 
and most of the H I mass, account for < 20% of the 
ionized mass in these clumps: due to the higher density 
of the gas in the tip and tail, the emission is relatively 
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bright, but L H + and, therefore, the H + mass are rela- 
tively small. 

The choice of assuming that the neutral and ionized 
gas are of equal pressure and temperature or equal den- 
sity results in a factor of two change in the derived mass 
but a factor of four change in the derived path length of 
ionized gas. For the tip and tail, the data (Table [2| are 
consistent with either equal pressure or equal density: 
the masses are small and the path lengths much smaller 
than our projected resolution in both cases. However, 
for both clumps A and B, equal pressure yields a path 
length (Xh+) considerably larger than the projected size 
of the clump on the sky, whereas equal density yields a 
value of i H + similar to the projected size of the clump. 
Therefore, unless the clumps are elongated along the line 
of sight by a factor of 2 — 4, an ionized skin in pressure 
equilibrium with a neutral clump does not effectively de- 
scribe clumps A or B. If the temperature of the neutral 
gas is lower than the temperature of the ionized gas, the 
pressure equilibrium ionized skin describes clumps A and 
B even more poorly. Our data do not rule out an ionized 
skin of the same density as the neutral gas in clumps A 
and B. Pressure effects should act in a sound crossing 
time, < 40 Myr for these clumps (assuming T = 10 4 K), 
although external dynamical forces could prevent a static 
equilibrium from being established in these more tenu- 
ous region s (Benjamin &; Danly 1997 Peek et al. 2007 
Heitsch fc Putman||2009 = | ~ 



definition. 



The dominant source of uncertainty within this model 
is the distance to the cloud. Because nj 1 oc ijj j cx D, 
the mass scales as M oc EMD 3 . The range of al- 
lowed distances quoted by L08 yields a ±30% statis- 
tical uncertainty in our mass estimate. Applying an 
extinction correction would increas e th e mass. If the 
A(H.a) = 0.38 mag we estimated in § 2.3 applies through- 
out the cloud, the mass estimate would be roughly 40% 
higher than stated here. 

Our derived ionized gas mass estimates are sensitive to 
the chosen boundaries of the neutral clump and the size 
of the clump, Ljj j. Mass estimates scale linearly with 
Ljj j, which we choose by eye based on the 21 cm emis- 
sion. We identify the clumps themselves subjectively, 
choosing regions with enhanced Ha emission and roughly 
corresponding enhanced 21 cm emission. For the tip, 
clump A, and clump B, this procedure is reasonably ro- 
bust; the H I contours are steep enough so that varying 
the chosen Min TVjj j by a factor of 2 changes the de- 
rived H I mass by ~ 40%. The 21 cm emission defining 
the tail is considerably less robust, as a number of sim- 
ilar enhancements in H I are evident in the body of the 
Smith Cloud (see Fig. |2|. We identify the cloud based 
on the strong Ha emission and choose the H I contour 
level so that the region with the neutral emission roughly 
matches the region with the ionized emission. However, 
there is relatively little H + mass in the tail, so the cho- 
sen contour level is not crucial to the total mass estimate. 
Alternatively, one could define the H I boundary of the 
clump as the extent of the WHAM beams included in the 
clump. This approach modifies the derived H + mass es- 
timates by < 30% in the ionized skin model and < 10% 
in the fully mixed model. Because the H I data have 
considerably higher angular resolution than the WHAM 
Ha data, we do not adopt this alternative H I clump 



3.1.2. Fully mixed ions and neutrals 



We now estimate the mass of the clumps assuming 
that the H + and H I in each clump are fully mixed or, 
cquivalently, L H + = Ljj j. This scenario applies if the 
clumps consist of a number of smaller, unresolved clouds 
of uniform den sity and the source of ionization penetrates 
or permeates ( jBland-Hawthorn et al. 20071 the entire 
clump. We determine the mean emission measure from 
the clump, (EM), and derive the root mean square elec- 
tron density in the clump: n e — ((EM)L^ 4 ) 4 / 2 . The ion- 
ized mass in the clump is lAmnn e L H +flD 2 . The mass 
of each clump in this model is shown in Table [2] with 
a total of w 3 x 10 5 Af Q of ionized gas in the identified 

clumps. This estimate scales as Ma EM 1 / 2 D 25 , so the 
uncertainty due to the distance is ±30% and a 0.38 mag 
extinction correction would increase the mass by 20%. 

3.1.3. Total mass 

With both the fully mixed model and the ionized skin 
model, we find that the two clumps with low H I col- 
umn density (clumps A and B) each have H + masses 
of 3 — 8 x 10 4 M Q . In the two clumps with higher 
H I column density (the tip and tail), the ionized skin 
model places much less mass in the tip and tail (each 
0.6-0.8 xl0 4 M©), whereas we find M H + « 4-8xlO 4 M 
in each of the four clumps in the fully mixed model. 

The clumps identified in our mass estimate do not ac- 
count for all the Ha or H I emission in the cloud. Using 
the techniques described above, we calculated the mass 
in each of the beams not associated with a clump as- 
suming (1) n = 0.01 cm~ 3 , with L H + derived from the 
Ha and n e = no or n e = no/2, as in the ionized skin 
model, and (2) i H + = 1 kpc, deriving n e from the Ha 
emission measure, as in the fully mixed model. We es- 
timated these values based upon the extent of the H I 
emission and typical H I column densities in the regions 
outside our identified clumps. In the fully mixed model 
(L H + = Ljj j) and equal density ionized skin (n e = n ) 
models, the total ionized mass in the diffuse sightlines is 
«3x lO 6 Af , considerably more than that present in the 
clumps. If the temperature and pressure of the ionized 
and neutral gas are equal, the resulting mass is a factor 
of 2 higher. However, we consider this model implausi- 
ble for this case because the path length along the line 
of sight, i H + = 4 kpc, is much larger than the observed 
transverse size of the Smith Cloud. These results are all 
sumarized in Table [21 

3.2. Temperature 

Line widths of Ha and the relatively narrow [S II] line 
allow us to constrain the temperature of the ionized gas, 
assuming the ion ized hydrogen and sulfur are fully mixed 
( |Reynolds|[T985| : 



T 



10 4 K 



VF H 



21.1km s" 1 



21.0 km s" 1 



- 0.072, 



(2) 

where Wh and Wg are the full widths at half maximum 
of the Ha and [S II] lines. This solution requires a distri- 
bution of nonthermal velocities with the most probable 



G 



TABLE 2 

Mass estimates for individual clumps 



Clump 



Min ATjt f 
(cm -2 ) 



(EM) 
(pc cm -6 ) 



WHAM 
beams 



Skin (n e = no/2) Skin (n e = no) 



Mixed (L 



(1O<*M ) 



L HI 

(Pc) 



(no) 
(cm -3 ) 



H+ 



(PC) 



M H+ b 
(10 6 M Q ) 



L H+ 
(PC) 



M H+ b 
(10 6 Af Q ) 



(n e ) 
(cm -3 ) 



M H+ ~ 
(10 6 M Q ) 



Tip 200 xlO 18 

Tail 90xl0 18 

Clump A 20xl0 18 

Clump B 20xl0 18 
"Diffuse 3 



0.88+0.05 
0.31±0.07 
0.51±0.15 
0.39+0.08 
0.11+0.01 



3 0.6 

6 0.2 

3 0.07 

6 0.19 

244 5. 



200 0.4 

100 0.4 

200 0.06 

300 0.04 

1000 0.01 



18 0.015 

8 0.012 

500 0.06 

900 0.14 

4000 7. 



5 0.008 

2 0.006 

140 0.03 

200 0.07 

1100 3. 



0.07 
0.06 
0.05 
0.04 
0.009 



0.05 
0.04 
0.04 
0.08 
3. 



a H I column densities are integra ted fro m +70 to +145 km s _1 . b Mass estimates in the ionized skin model scale as EM and, therefore, linearly 

: Mass estimates in the fully mixed model scale as EM 1 / 2 and, therefore, as the square root of any 
refers to sightlines not included in the tip, tail, clump A, or clump B. 



with any extinction co rrectio n (< 
extinction correction (§ 3.1.2 I. ' 



3.1.11 



Diffuse" 




Fig. 4.— Nitrogen abundance ([N/H] = log(N/H) - log(N/H) Q ) 

derived from the observed line ratio, [N II]/Ha = 0.32 ± 0.05, as 

a function of the electron temperature, assuming N+/N = H+/H. 

The error envelope due to the uncertainty in the line ratio is shaded 

grey. The dot-dashed curve shows the nitrogen abundance derived 

from the P03 observation of [N II]/Hq = 0.6. Solid lines show the 

range of temperatures allowed by the line width analysis, while 

dashed lines show the range of typical temperatures in a diffuse 

gas heated by photoionization. The dotted line shows the solar 

abundance, 
value 



km s 



1.64 km s" 



! km s 



(3) 



The observed widths (Table 1) yield T = 15500 ± 7400 K 
and v = 11.5±3.4km s _1 in the direction of the long ex- 
posures in the tip. The large uncertainty is primarily due 
to blending of the the Smith Cloud and Sagittarius Arm 
emission. If the Sagittarius emission is better physically 
represented by multiple components, our three-Gaussian 
fit overestimates Wh in the Smith Cloud. With a smaller 
Wh alone, the temperature of the cloud would be lower. 
The red wing of the Smith Cloud Ha spectrum has nearly 
zero contribution from the Sagittarius Arm and there- 
fore provides a useful check on the width of the Smith 
Cloud component. Keeping the Ha mean fixed to the 
+96.2 km s^ 1 [S II] mean, line widths Wh < 27 km s _1 
fail to reproduce this portion of the spectrum. This yields 
a lower bound of T = 7000±2000 K, with the uncertainty 
propogated from the uncertainty in Ws- 

3.3. Nitrogen abundance 

We now use the pointed [N II] and Ha observations 
to constrain the metallicity of the Smith Cloud. The 
line ratio [N II] A6583 /Ha is sensitive primaril y to the 
temperatur e and nitrogen abundance of the gas ( Haffner 
eTaL][l999] |: 

H+\ _1 /N+ N 

IT/ I InTh 



[N H] 
Ha 



= 1.63 x 10 s 



T 



0.426 -2.I8/T4 



Because of the similar first ionization potentials of nitro- 
gen (14.5 eV) and hydrogen (13.6 eV), N+/N « H+/H; 
photoionization mode lling also supports this argument 
i |Sembach et"aL||2000[ ) . The observed [O III] upper limit 
of 0.07 K (ionization potential of 35 eV) supports an 
assumption that there is little N ++ (ionization poten- 
tial of 54 eV) in this gas. Therefore, the observed 
line ratio [N II]/Ha = 0.32 ± 0.05 constrains the nitro- 
gen abundance, shown in Figure [4] The lowest allowed 
temperature, s» 8000 K, yields a nitrogen abundance of 
N/H = 0.4 times s olar (assuming (N/H) Q = 7.5 x 10~ 5 ; 
Meyer et al.|1997| , whereas the highest allowed temper- 
ature, « 23000 K, yields N/H = 0.05 times solar. The 
figure also shows the range of typical temperatures in 
the WIM, which better constrain the nitrogen abundance 
with assumptions discussed below. 

4. DISCUSSION AND FUTURE WORK 



Typical Ha intensities for HVCs are 0.06-0.5 R flTufte 
et al.||1998||Haffner et al 120011 [Werner et al.|2001[ 



Putman et al. 



(4) 



2003| |Haffner||2005P ; the intensities ob- 
served here for the Smith Cloud are in this range. The 
narrow range of Ha intensities for these clouds suggests 
that, like the WIM, all of these HVCs are ionized by a 
diffuse UV radiatio n field, not shocks or individu al, local 
hot stars (see also |Bland-Hawthorn et alT 19981, which 
would produce significant variations in Ha intensity from 
cloud to cloud and even within a cloud. The similar line 
ratios [N II] /Ha = 0.32±0.05 and [S II] /Ha = 0.32±0.05 
observed here are co mmon in WIM gas b ut atypical of 
classical H II regions dMadsen et a^[2006l ). 

The extent of the Ha emission from the Smith Cloud 
generally traces the H I, but we find moderately bright 
Ha emission in regions with faint H I. This result is sim- 
ilar to that for the intermediate velocity clouds Com- 
plex K and Comple x L, which have also been mapped 
in Ha with WHAM ( |Haffner et aTpjOll |Haffner|[2005| . 
In all three clouds, the presence of Ha and H I emission 
tend to track each other in corresponding sightlines, but 
there is no strong correlation between the Ha and H I 
intensities, with the exception of the tip of the Smith 
Cloud (blue ' x ' in Fig. |2| , which is the brightest portion 
in both Ha and H I. Note that the peak Smith Cloud H I 
column density of 5.2 x 10 20 cm -2 (L08) is a fact or of 10 
higher than that for Complex K, 5.3 x 10 19 cm -2 ( |Haffner| 

" (P03). The 



|et al.||200l| , or Complex L, 3.6 x 10 
lack of a strong correlation between Ha intensity and H I 
column density in the low-column density region of the 
cloud supports previous suggestions that the ionized gas 



if) 



cm 
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is not mixed with the neutral gas and that variations in 
the Ha intensity outside the high-H I column density tip 
may be due to variations in the strength of the ionizing 
radiation field rather than the distribution of neutral gas. 
The Ha maps of Complex K, Complex L, and the Smith 
Cloud all suggest that ionized gas extends beyond the 
H I associated with the clouds. Moreover, there is some 
evidence for high velocity gas with higher ions, including 
C IV, along a sightline with no ev idence of lower ions or 
neutral gas (Sembach et al.| |1999[) . 



Our analysis of the Ha and |S II] line widths leaves 
th e te mperature of the Smith Cloud very uncertain 
(§ 3.2 1. We estimate constraints on the temperature 



by comparing the physical conditions in the HVC and 
the WIM. Typical observed temperatures of WIM gas 
are « 6000 — 10000 K, with h igher temperatures in lower 



density por tions of the WIM (Haffner et al. 1999 Madsen 
et al.||2006|). Models show that photoionization heating 



from diffuse UV radiation balanced by radiative cool- 
ing does not account for the observed temperatures, sug- 
gesting that an addition al heating source is important 



at densities < 0.1cm 3 (Reynolds et al. 1999). Also, a 
low metallicity would reduce the dominant cooling mech 



anism, collisionally excited line emission (Oster brock & 



Ferland 2006); the WIM has roughly solar metallicity. 
The low density and low metallicity of the Smith Cloud 
suggest that the temperature of the ionized gas in the 
cloud is > 8000 K, the temperature of relatively warm 
WIM gas. The implied nitrogen abundance is then 
0.44 ± 0.07 or 0.15 ± 0.02 times solar assuming T = 8000 
or 12000 K, respectively (see 



:-!.:-! ,nul Fig.pJ. 



For com- 
0.2 times 



pariso n, the metallicit y of Complex C is 
solar ( Fox et al. 2004 and references therein) . 

P 03 found a line ratio [N II] A6583/Ha = 0.60 in each 
of their fields, significantly higher than the value we find 
in the 1° WHAM beam in the tip of the cloud. We spec- 
ulate on three possible explanations for this difference: 
(1) A discrepancy in the data taking, reduction, or anal- 
ysis, (2) effects of the different beam diameters, 1° in 
our WHAM data, compared to 5' in their data, or (3) 
variations in the nitrogen abundance across the cloud, 
as the P03 data were taken approximately 2° down the 
major axis of the cloud from the tip, where our pointed 
spectra were obtained. Beam effects are possible if there 
is clumping on angular scales not resolved by WHAM. 
With the typical assumption that individual HVCs are 
coherent objects with uniform metallicities, a variation 
in the line ratio along the axis of the cloud would mean 
that the tip of the cloud is cooler than the downstream 
portion observed by P03. 

Our mass calculations (§ |3.1[ ) in combination with a 
comparison of the Ha and HI morphology of the cloud 
suggest that a different picture is appropriate in the 
area of the cloud along the major axis (including the 
tip and tail) than in the area of fainter emission (includ- 
ing clumps A and B). Along the major axis, the data are 
consistent with an H + skin on an optically thick (to ion- 
izing photons) H I cloud: the derived skin depth is small 
(~ 10 pc) and the Ha emission does not extend consid- 
erably beyond the H I. Meanwhile, in clumps A and B 



and assuming that the H + thermal pressure equals the 
H I pressure, the required thickness of the ionized layer 
is considerably larger than the extent of the observed 
Ha emission, implying that the ionized and neutral gas 
are mixed. This picture indicates that there is relatively 
little ionized gas associated with the brightest Ha and 
21 cm emission. The total H + mass in the region we 
have observed is ~ 3 x 1O 6 M (Table The distribu- 
tions of the H I and H + mass are different: over 80% of 
the H I is along the major axis of the cloud, above the 
green line shown in Figures [2] and [3j while the H + masses 
above and below the green line are approximaely equal. 
Because Ha emission continues to the edge of our map, 
we cannot say how far the Ha emission (and H + ) extends 
beyond the 21 cm emission. 

The kinematics and cometary morphology indicate 
that the Smith Cloud is moving towards the Milky Way. 
The immediate question is whether the cloud traces new 
infall or gas cycling back to the Galaxy as part of a 
"Galactic fountain". L08 argued that, given the trajec- 
tory and velocity of the cloud, it is unlikely that the 
material was ejected from the disk, although the low in- 
clination and prograde orbit of the cloud suggest that the 
orbit has been affected by the Galactic potential. Our 
metallicity constraint supports a non-fountain origin. 

WHAM began operations from Cerro Tololo Inter- 
American Observatory in Chile in April 2009. As the 
Smith Cloud is visible from this southern site as well, 
we will obtain follow-up observations of the region to ex- 
pand upon the results presented here. Specifically, we 
plan to: (1) Extend the Ha map beyond the limits pre- 
sented here. Ha emission from the Smith Cloud extends 
to the edge of the region we have currently mapped. (2) 
Obtain H/3 observ ation s that allow more accurate extinc- 
tion correction (§ 2.3), which will improve the estimate 



of the ionized mass m the cloud. (3) Improve the signal- 
to-noise ratio of the [O III] observations to obtain more 
stringent constraints on the ionization conditions within 
the cloud. (4) Add additional, deep observations of Ha, 
[N II] , [S II] , and [O III] at multiple points in the cloud 
to improve our understanding of the kinematics of the 
ionized component and the physical conditions through- 
out the cloud. Tracing variations in these values should 
provide considerable insight into the interaction of in- 
falling material with the disk of the Galaxy. In particu- 
lar, multiple observations of the relatively narrow [S II] 
line may allow us to identify variations in the velocity 
of the ionized gas; we could not identify any variations 
in the velocity of the S mith Cloud emission with the Ha 
maps presented here (§ 2.2 1. 
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